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' p Abstract 

. We consider the ortho- and parapositronium production in the process AA 



AA+ Ps where A is a nucleus with the charge number Z. The inclusive cross 
section and the energy distribution of the relativistic Ps are calculated which 
' are of primary interest from the experimental point of view. The accuracy of 
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the corresponding cross sections is given by omitting terms ~ (Za) 2 /L 2 for 
the para-Ps and ~ (Za) 2 /L for the ortho-Ps production where L = Iwy 2 ~ 9 
and 16 for the RHIC and the LHC. Within this accuracy the multiphoton 
(Coulomb) corrections are taken into account. We show that the RHIC and 
the LHC will be Ps factories with a productions rate of about 10 5 -j- 10 8 
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relativistic Ps per day. The fraction of the ortho-Ps is expected to be of the 
same order as that of the para-Ps for Au-Au and Pb-Pb collisions. 

PACS numbers: 25.75-q, 36.10.Dr, 12.20.-m 
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I. INTRODUCTION. 



Parapositronium (para-Ps, n 1 Sq) has a positive charge parity C = +1 and in the ground 
state, n — 1, its life time at rest To is cro = 3.7 cm. On the other hand, orthopositronium 
(ortho-Ps, n 3 S\) has C — — 1 and cro = 42 m. The relativistic Ps has a decay length 
CToe/rrips where e and mp s « 2m e denote its energy and mass. The production of Ps in 
different collisions is of interest by the following reasons: 

• Ps can be used to test fundamental laws like the CPT-theorem (see review [|T|). 

• A more detailed comparison between the experimental value of ortho-Ps width and 
its theoretical prediction is necessary since up to now there is an essential difference 
between them 0). 

• Relativistic Ps can be used to study the unusual large transparency of Ps in thin layers 
predicted in Refs. |||. 

Up to now the following methods have been proposed to create positronium: The photo- 
||- |§ and electroproduction |5|,|9|Jirj|| of Ps on nuclei and atoms. The results of Refs. ||- fi~0[ 



were completed and corrected in Ref. [11]. Recently a new method, based on an electron 
cooling technique, was proposed in the works 0,0. It was claimed that this method can 
give about 10 4 nonrelativistic Ps per second on a specifically designed accelerator. 

In the present paper we study in detail a new and promising possibility — the production 
of relativistic Ps at colliders with relativistic heavy nuclei like the RHIC and the LHC, i.e. 
we consider the process AA —>■ AA + Ps where A denotes the nucleus with charge number 
Z. The parameters of these colliders (from Refs. [13]) are given in Table |. We show that 
at these colliders 

~ 1 -T- 10 3 relativistic Ps per second (1-1) 

will be produced.lt is interesting to note that the ratio of ortho-Ps to para-Ps events turns 
out to be large: for the RHIC about 0.8, for the LHC 0.3 in the Pb-Pb mode and 0.02 in 
the Ca-Ca mode. 

Up to now only the total cross section for para-Ps production has been calculated in the 
leading logarithmic approximation (LLA) (see Ref. ||). There it has been found: 

< L r t = 7rC(3)a L 3 , a = — ^, ((3)^4 = 1.202..., L = ln 7 2 . (1.2) 

m e n=l U 

However, for experiments (especially with relativistic Ps), it is much more interesting to know 
the energy and angular distributions of Ps which are given by the inclusive cross section and 
energy spectrum of Ps. These distributions for para-Ps as well as for ortho-Ps are calculated 
in the present work. In our calculations we also take into account with a certain accuracy the 
Coulomb correction (CC) related to multiphoton exchange of the produced Ps with nuclei 
described by Feynman diagrams of Fig. |l[ They are quite important and decrease the cross 
sections by 30 4- 50 % for the Au-Au and Pb-Pb collisions. The approach used here is based 
essentially on results on the production of relativistic positronium in collisions of photons 
and electrons with heavy ions presented in our recent work JTT . 
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Our main notations are given in Eq. ( |1.2|) and Fig. ^ the i-th nucleus changes its 4- 
momentum (energy) Pj (E) to P( (E'j) and emits n or n' virtual photons with the total 
4-momentum q^ = Pi — P[ and the energy cu, = E — E[. These photons collide and produce 
Ps with 4-momentum p = qi + q 2 and energy e = ui\ + uj 2 - Due to charge parity conservation, 
the sum n + n' is even (odd) for para-Ps (ortho-Ps) production. We use the c.m.s. with the 
z-axis along the vector Pi and denote by p_i_ the transverse components of vector p. The 
CC depends on the parameter 

v = Za^ Z/137 . (1.3) 

In Section 2 we discuss the production of para-Ps in ultrarelativistic heavy ion collisions. 
After selecting the leading diagrams including the multiphoton exchange with the nuclei, 
the Born cross sections are calculated in the equivalent photon approximation to identify the 
dominant kinematical regions. The exact Born cross sections for ultrarelativistic ions are 
found in the next subsection. The CC are calculated in the accuracy given by the selection 
of the leading diagrams. Section 3 is devoted to the production of ortho-Ps. Our results are 
summarized in Section 4 where we discuss the expected production rates for the colliders 
RHIC and LHC. 



II. PRODUCTION OF PARA-PS 

A. Selection of the leading diagrams including the multiphoton exchange 

The amplitude M para of the para-Ps production can be presented in the following form 
via the sums over amplitudes M nn i of Fig. 1 (with n + n' being even) 

M para = Y M nn , = M B + M 1 + M l + M 2 , M B = Mh, (2.1) 

n,n'>l 

Mi = Yl M m> , Mi = M nl> M 2 = Y M nn> ■ 

n'=3,5,7,... n=3,5,7,... n,n'>2 

Here Mb is the Born amplitude in which each nucleus emits one photon, the item Mi (Mi) 
corresponds to the case when the first (second) nucleus emits one photon while the other 
nucleus emits at least three photons, the last term M 2 corresponds to the case when both 
nuclei emit at least two photons. Using this decomposition of the amplitude, we present the 
cross section in the form 

0"para = + 0"l + + ^2 , d<J B OC |M B | 2 , 

do-i oc 2Re (M B M*) + |Mi| 2 , da 1 oc 2Re (M B M*) + |M| 2 , (2.2) 

da 2 oc 2Re (M B M* + M X M[ + M X M 2 * + MM) + |M 2 | 2 . 

It is not difficult to estimate that the integration over the total virtual photon momenta 
squared q\ and q\ gives two large Weizsacker-Williams (WW) logarithms ~ L 2 in ctb, one 
large WW logarithm ~ L in o\ and &\ and no one large WW logarithm in cr 2 . Therefore, 
the relative contributions of the items <Tj is 0\j a B = a\j a B ^ u 2 /L and 
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2 

2l „ 2- < 0.4% (2.3) 

for the colliders in Table |. As a result, with an accuracy of the order of 1 % we can neglect 
cr 2 and use the equation 

a = a B + 0i + crj . (2.4) 

Moreover, with the same accuracy we can calculate o"i and d\ in LLA only since the terms 
next to the leading ones are of the order of v 2 jL 2 . 

B. Calculating the Born cross section in the equivalent photon approximation 

Before we proceed to a more exact calculation of da B , it is useful to perform a rough 
estimate using the equivalent photon approximation (EPA). This allows us to identify the 
dominant region in phase space. 

It is well known that the main contribution to the pair production AA — > AAe + e~ is 
given by the region of small photon virtualities — q 2 <C (qi + q 2 ) 2 - Therefore, for the Ps 
production this region isQ 

- Qi « 4m 2 , (2.5) 
i.e. the virtual photons are almost real in the dominant region. Since 

-% 2 = q 2 ± + (^/7) 2 , (2.6) 

we immediately conclude that in the dominant region the photon energies u>i are small 
compared with the nucleus energies E: 

Ui ~ m e 7 < E . (2.7) 

In the EPA the cross section for the AA — > AA + Ps process is related to the cross section 
for the real photoprocess 77 — ► Ps by 

da B = dn^ujt) dn 2 {uJ2) 077-^(^1^2) , u 1)2 ~ -{e±p z ) , (2.8) 

where drii is the number of the equivalent photons generated by the i-th nucleus. For drii 
we use two equivalent forms (see, for example, formulae from Appendix D of review [|14j] in 
which we omit terms of the order of u>i/E (in virtue of relation (|2"7! 

Z 2 a q 2 ± du, 
7T 2 (q 2 ) 2 UJ, 



dn l = ^ r j^ — d 2 q l± , (2.9) 



This can also be seen from the expression for the Born amplitude obtained in Appendix A. 
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dm 



Z 2 a 



71 



1 - 



du}j d(-qf) dpi 
Ui (-q 2 ) 2tt 



(2.10) 



where ifii is the azimuthal angle of the vector q,. It is also well known that the equivalent 
photons of small energy are completely linearly polarized in the scattering plane and their 
polarization vectors are 



(2-11) 



The cross section for the real photoprocess 77 — > Ps can be easily obtained using the 
known width of the n 1 Sq state Tp s ^ 77 = a 5 m e (ei x e 2 ) 2 / {An 3 ) and summing over all para-Ps 
states n 



(T 77 ^ Ps = 4vr 2 a 5 C(3)(e 1 x e 2 ) 2 5(Auj l uj 2 - Am 2 ) . 



(2.12) 



Substituting Eqs. (|2.10|) - (|2.12| ) into ( |2.8|) and integrating over <pi, oj 2 and over qf in the 
region 



2 < 2 < 2 



(2-13) 



we find with logarithmic accuracy 

da-B = 2(, (3)cr In In , lo\UJ 2 ~ m e . (2-14) 

Taking into account the inequalities (|2.7|) , we conclude that the main region for u\ is 

— ~ c<Ji ~ m e 7 . (2.15) 

Integrating rfo"B ( P-14 ) over the photon energy u)\ in the region (|2.l5j ) we obtain the total 
cross section (|1.2j) . 

The spectrum of the relativistic Ps can be obtained from Eq. ( |2-14 ) as follows. Let us 
divide the region ( 2.15|) into two symmetrical subregions. In the first one, m e <C iO\ ~ m e 7, 
we have uj\ ~ e ~ p z , uj 2 ~ m 2 /s and the contribution to the spectrum is 



d4 1) = 2C(3)a ln^ln£l^ 



(2-16) 



In the second subregion m e /j ~ u>i <C m e we have u 2 ~ £ ~ — p z , t^i ~ ^e/ £ an d 



2 C (3)a ln£lln^^. 

m e £ £ 



The total contribution of these two subregions is 



a A/-fo\ 1 m ^ 1 de s~ < 

do b = 4(73) (To In In , m e <£~ w e 7 • 



(2.17) 



(2.18) 



Integrating this spectrum over £ leads to the result ([TTJ 
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To calculate the inclusive cross section we substitute Eqs. (|2.9|) , ( |2.11|) , ( |2.12| ) into 
(2.8) and perform a simple integration 



J uji u 2 4m2 e 

Further, we insert into the right-hand side of Eq. ( |2.8| ) the factor 

1 = J S(q 1± + q 2 ± - pi) d 2 p ± 

and find the inclusive cross section in the form 

d 3 c?B C(3) f (qu x q 2i ) 2 2 2 

^ "d^" = ^ J (g g 2)2 <^ + ^ - P-L) d <^ ?2± • ( 2 - 19 ) 

From this expression and taking into account Eqs. ( |2.6|) , ( |2.13|) , we conclude that the leading 
logarithmic contribution to the cross section is given by the region 

/ \ 2 / \ 2 

/ Wi \ / UJ2\ 2 ^ 2 



V 7 ) \ 7 . 
It arises from the two symmetrical subregions: 



<pl<m 2 . (2.20) 



<q?±<pi' q2±~p±, 

— 1 <qL<PL qi±«P±- (2.21) 



The contribution from the first region is 



C(3) ^ , Pi , 2 (qi± x P± ) 2 C(3) cro p 2 ±7 2 



2tt 2 7o (q?±Pl) 2 2tt p^ " w 2 



Taking into account the second region, we find the whole inclusive cross section 

d 3 a B _ C(3) a ln pj 7 2 
G? 3 p 7r pj_ m 



e^ = ^^ln^f. (2.22) 



Integrating this expression over p^ in the range ( 2.20| ) and taking into account the two 



symmetrical subregions ( |2.16| , |2TT7| ) mentioned before, we obtain the spectrum ( |2.18| ). 



C. The exact calculation of the Born cross section for ultrarelativistic ions 

The Born amplitude of the discussed process is obtained in Appendix A omitting terms 
of the relative order of I/7 2 , m e /M or smaller. Using this expression for the amplitude we 
can present the inclusive cross section in the form 

e = C(3) ^± j B ; (2.23) 
d 3 p An m 2 ' 
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Am? f 

J B = - / A 2 5(qi_L + q 21 _ - p ± ) d 2 q 1± d 2 q 2 ± , 

71 J 

A = qi± x q 2± 4m 2 

q\ q\ Am\ ~ Qi ~ Q2 

which differs from the approximate Eq. ( |2.19| ) only by the factor (Ami) 2 / (4ml ~ Qi ~ Q2) 2 
in the integrand. It is convenient to introduce the dimensionless variables 

x = - , x = , t—[ . (2.24) 

4m e 7 4m e 7 \2m e J 

They are not independent, but connected by the constraint 

1 + r , 

In this notation 

-ql = q 2 l± + (2m e x) 2 , -q\ = q 2 ± + (2m e x) 2 . (2.26) 

It is useful to note that the rapidity of the Ps is equal to (1/2) ln(x/x) and that the 
variables x and x have a very simple form for relativistic Ps. Indeed, for p z 3> m e the variable 
x is x = e/(2m e , y), i.e. x is the ratio of Ps energy to the characteristic energy 2m e r ) (above 
2m e 7 the spectrum drops very quickly). Analogously, for (— p z ) ^> m e x = e/(2m e r y). The 
neglected contributions are of the order of (4m 2 + pl^/e 2 . From the discussion in Sect. [II B 



it follows that the dominant contribution to the cross i section ( [2.23| ) is given by the region 

— ~ x , x ~ 1 , x 2 + x 2 ~r~l. (2.27) 

T 

In Appendix B we show that the function J B can be reduced to the one-dimensional integral 

r r 1 

J B = - u 3 (l-u)(2-u) 2 K(u,x,T)du, (2.28) 
8 Jo 

f 1 (l-t 2 )dt 

Ar(u,x,r)=3/_ i ( J 2 + M ; c)4 (2-29) 

12a(o 2 - 5a 2 + ac) , c - a + \fd AN 
In -= + 



where 



dV 2 c-a-Vd [(a + c) 2 -6 2 ] 2 rf 3 

N = -30a 6 + 6 6 - 104a 5 c + 56a 4 6 2 - 124a 4 c 2 + 90a 3 6 2 c - 56a 3 c 3 
-26a 2 6 4 + 36a 2 6 2 c 2 - 6a 2 c 4 - 2a6 4 c - 6a6 2 c 3 



a = —-u r , b = - u(2 — u) (x — x ) , 
c= -(2-u) 2 t + \(2-u) 2 (x 2 + x 2 ) + (l-u)(2-u), d = b 2 -Aac. (2.30) 



S 



Let us stress that the inclusive cross section ( 2.23|) with ( 2.28|) is valid both for relativistic 



para-Ps with e 3> 2m e as well and for nonrelativistic para-Ps with e ~ 2m e . 

Using the representation ( |2.23| ) for J B it is not difficult to find that the inclusive cross 
section is large in the main region 



T 



ln(4 7 2 r) - 1 at x 2 + x 2 < r < 1 (2.31) 



but drops very quickly outside this region. Indeed, the function Jb is small at large r 



J B = \ ln(4 7 2 ) at x 2 + x 2 < K r , (2.32) 



at large a; 



J B = -7 ln(4 7 2 a; 2 ) at r < 1 « x 2 , (2.33) 

x 8 



and at very small r 

Jb = 7 2 ' - 2 . 2 I " ' " In ^ - 1 ] at r < x 2 + x 2 < 1 . (2.34) 




In Figs. we present the inclusive cross section ( |2.23| ) at the RHIC and the LHC 



(Pb-Pb mode) as function of the scaled Ps energy for different scaled transverse momenta 
squared r. From these Figures it is obvious that for e/{2m e ^) < ^fr the cross sections only 
weakly depend on the Ps energy. The cross sections quickly drop at larger energies. 

To obtain the spectrum of relativistic Ps it is convenient to use an approach developed in 
Ref. [IlJ for the electroproduction of Ps on nucleus eA — > eA + Ps (Fig. ^). It is well known 
that the electroproduction cross section can be exactly written in terms of two structure 
functions or two cross section 0"r(s 7 , q 2 ) and crs(s 7 , q 2 ) for the virtual processes j^A — > Ps+A 
and jgA — * Ps + A (where 7^ and 7J denote the transverse and scalar virtual photons with 
4-momentum q and helicity Xt = ±1 and A5 = 0, respectively, and s 7 = 2qP 2 ): 

da(eA — > Ps + A) = cry(s 7 , q 2 ) dnr(s 7 , q 2 ) + o"5(s 7 , q 2 ) dns{s 1 ^ q 2 ) . (2.35) 

Here the coefficients dnx and dns can be called the number of the transverse and scalar 
photons generated by the electron. 

To calculate da B in our case we can use a similar expression where the number of photons 
are now generated from a scalar projectile nucleus. i The cross sections <jt and as for para- 
Ps have been obtained in Ref. |IT] for high energetic virtual photons (for s 7 ^> 4m e M or 
7 2 x > 1): 

Z 2 a e L(x,y)-1 

a T = ttC(3) — — rrj-, 0-5 = (2.36) 



with 



ml (l + y)' 



L(x,y) = ln(4 7 2 x)-^ln(l + t/), y = ^ . (2.37) 
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Here x is defined in Eq. (|2.24j ) . The accuracy of Eq. Q2.36Q is determined by neglecting terms 

~ l/(7 2 a;), y—q%/(m e j 2 x). The number of photons generated by the first nucleus are (see, 
for example Appendix D from Ref. 



drti 



Z 2 a 



1 - 



x 2 \ dx dy 



dn<- 



it \ y j x y 

Z 2 a dx dy 
Ti x y 



Note that driT corresponds to dn\ as given in Eq. ( [2.10| ). 

Using Eqs. (|2.36|) -( [2~38D and integrating over y we obtain^ 



do\ 



C(3)a {hix) [ln(4 7 2 a;) - l] -/„(*)} 



dx 



x 



with 



dy 



00 I x*\ 

2 V y) v( l + y) 2 



/ ,\ , 1 + X 2 

(l + 2x 2 ) In — — n 2 . 



ln(l + y) dy 



l + 2x 2 



IT 2 . ( X 2 ' 

~6 ~ Ll2 I 1 + x 2 



1 2 1 1 + x<2 
x in — 



- ln(l + x 2 



Here the dilogarithmic function is defined as 

ln(l - t) 



Li 2 (x) 



dt. 



For small x < 1 we have 

h{x) = 

while for large x ^> 1 



2 In 



x 



Mo) 



7T 2 -6 

12 



0.3225 



1 

6V 



1 A 5 
— - IniH 

6x 4 V 12 



^2.38) 



(2.39) 



(2.40) 



(2.41) 



(2.42) 



^2.43) 



In obtaining Eq. ( |2.39 ), we have considered the first nucleus as projectile and the second 
nucleus as target. Taking the first nucleus as target and the second nucleus as projectile we 
find for da-Q the same expression as in Eq. ( |2.39|) with the replacement x — > x: 



da B = C(3) a (fax) [ln(4 7 2 £) - l] - / 2 (*)} ^ . 



[2.U) 



2 The lower limit is y m in = x ^ \ the upper limit can be extended to infinity due to the fast conver- 
gence of the integral. 
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To obtain the spectrum da B /de we have to take into account the following two circum- 
stances for the relativistic Ps: i (i) the variables x and x are related in a simple way to e; (ii) 
there are two symmetrical regions in the c.m.s.: p z ~ e ^> 2m e and (— p z ) ~ £ ^> 2m e . In 
the first region the variable i is i = e/ (2m e 7) and we have to use Eq. (|2.39|) , in the second 
region x = £/(2m e7 ) where Eq. ( j2.44j ) is used. As a result, the spectrum for the relativistic 
Ps has the form (compare with the approximate Eq. (|2.18 )) 



da B = 2C(3) cr 



h 



2m e 7 



In 



2£ 7 



I -/2 



2m e 7 



de 
e 



(2.45) 



The accuracy of this expression is determined by the omitted terms of the order of I/7 2 , 
2m e /(e 7 ), (2m e /e) 2 . 

Note that in the main region 2m e < e < 2m e7 the spectrum is logarithmically enhanced 



da B = 4C(3) cr 



ln^ln^-ln(4 7 2 ) 



30-7T 2 

24 



de 
e 



(2.46) 



while for large e 3> 2m e7 the spectrum drops very quickly 



da-t 



C(3) /2m e7 ^ 4 

T a ° It 



ln(4 7 2 ) - g 



(2.47) 



The spectrum of the para-Ps is shown in Fig. |5|. The solid and dash-dotted lines present 
the result for the RHIC and the LHC (Pb-Pb mode) obtained by a numerical integration of 
the inclusive cross section ( |2.23|) , ( |2.28|) . The dashed and dotted lines present the spectrum 
for both colliders calculated from Eq. (|2.45 ) which is valid for relativistic Ps with e ^> 2m e . 
It can be seen that in the region 10m e < e < (0.1 -j- 0.2) 2m e7 the quantity edo^jde has a 
very weak dependence on e whereas the spectra drop quickly for £/(2m e7 ) ~ 0.5. 

From the experimental point of view, it is interesting to integrate the cross section over 
the region e > e ni \ n (where 2m e £ min 2m e7 ). It can be obtained from Eq. ( |2.45| ) 
omitting pieces of the order of £ 2 nin /[(2m e7 ) 2 /]: 



cr B = 2C(3) (T 



ln(4 7 




21 + 



7T 



12 



12 



2, 30- 

-I s H 

3 12 



7T 



I 



(2.48) 



where 



7T 2 + 3C(3) + 21 
c = — = 2.873 . 



/ = In 



2m e7 



Using this expression for £ min = 15m e and the results of a numerical integration in the 
region 2m e < e < 15m e , we obtain the total cross section for the Born contribution: 



a B = 17.8 mb for RHIC , 110 mb for LHC (Pb-Pb) , 0.42 mb for LHC (Ca-Ca) . (2.49) 
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D. An approximate calculation of the Coulomb correction 



Let us consider the correction do\ defined in Eq. ( |2.2|) . Its contribution to the spec- 
trum can be calculated just by the same method as it was done for the Born contribution. 
According to Ref. [|ll|] we have to replace in Eq. (|2.36|) 



L(x,y) 



-C(iz) 



from which it follows 



d<ri = -C(3) a C{v)h{x) 



dx 



x 



The quantity C(u) is defined with the help of function F(a, b; c; z) 
metric function: 



(2.50) 
(2.51) 

the Gauss hypergeo- 



C(iz) 



F(r,u) dr , z 



F(t, u) 



2r 1 + r 



1 + r 

7TV 



l2 



F(iv,-iv;l;z, . 

sinn txv 



^2.52) 



The curve C(y) and properties of the function F(r,u) are given in [|Tl|], we mention here 
only that C(u) = 0.6734 for Au, 0.7080 for Pb and 



Analogously, 



C{u) = 2.9798 v 2 at v 2 < 1 . 



<& 1 = -C(3)(ToC(y)Mx) — 

x 



(2.53) 



(2.54) 



As we have explained in Sect. |II A| , the spectra da\/de and d&i/de have to be calculated 
with logarithmic accuracy only. Within this accuracy we have x ~ w 1 /(2m e 7) and x ~ 
6o>2/(2m e 7). Then for do\ we can repeat the calculations similar to those which we have 
already performed to get the spectrum ( |2.18| ) from expression ( [2.14| ). In our case this gives 



da 1 = -2C(3) cr C(y) m 7 ; 



ds 



^2.55) 



Taking into account the same contribution from du\ we find (with logarithmic accuracy) 



de 

d{a Y + <7i) = -4C(3) cr hi7 2 — , 2m e < £ < 2m e 7 . 
The total contribution of CC to the para-Ps production is 



a 



cc 

para 



a x + ox = -2C(3) a C(z/)[ln 7 



2i2 



from which it follow that the relative contribution of CC is 



a 



CC 
para 



LLA 

para 



O 



6C(v) 
In j 2 

12 



(2.56) 



(2.57) 



^2.58) 



i.e. it is -43 % for RHIC, -27 % for LHC (Pb-Pb) and -2.2 % for LHC (Ca-Ca). 

To obtain the contribution of do\ to the inclusive cross section we can explore Eq. ( p. 52 ). 
Strictly speaking, the variable r in that equation is r = p^/(4mg — qf). But in LLA used 
here for da\ we can take the function F(r, v) on the mass shell of the first virtual photon 
(at q\ = 0), i.e. we can use r = pj_/(4mg). Now with the help of a simple relation 



dx 



x 



dr 



d 3 p 



Anm? e e 



we obtain 



ggi _ C(3) (To 
d 3 p 



An m\ 



J u J x = -h{x)F(T,u) 



(2.59) 



(2.60) 



Analogously, for d&\ we find 



Ji = -fi(x)F(r,v). (2.61) 
The function fi(x) is defined in (|2.40|). As a result, the inclusive cross section reads 



para 



d 3 p 



03) ^ T T 

i n •'para ; '-'para 

47r mi 



JB-\h(x) + h(x)}F(r,u) 



(2.62) 



where Jb is defined in ( |2.28| ), x and x are defined in Eqs. ( |2.24j) . In the LLA we have 

fi(x) + fi(x) = 21n 7 2 , 2m e <£<2m e7 . (2.63) 



III. THE PRODUCTION OF ORTHO PS 

The amplitude of the ortho-Ps production can be presented in the form (compare with 
Eq. O) ' 

M ortho = M 1 + M 1 + M 2l (3.1) 
M 1= E M m>, Mi= E M -l' M 2= E M nn'- 

n'=2,4,6,... n=2,4,6,... n,n'>2 

with sum n + n' being odd. The total cross section reads (compare with Eq. (|2.2|) ) 

o"ortho = o\ + o\ + cr 2 , d(7i oc |Mi| 2 , d5i oc \Mi\ 2 , (3.2) 

rfa 2 oc 2Re (M 1 M* + + M X M^) + \M 2 \ 2 . 

Since Re (M\M*) disappears after azimuthal averaging, it is not difficult to estimate that 

- ~ < 4% ( 3 - 3 ) 
for the colliders in Table [I]. Therefore, 
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Cortho = cr 1 + a\ (3.4) 

with an accuracy of the order of 4 %; moreover, with the same accuracy we can calculate 
o\ and ax in LLA only. 

To obtain the inclusive cross section we use the relation (|2.35|) with driT and dns from 
Eqs. ( |2.38| ) and the cross sections d<jT and da$ for ortho-Ps from Eqs. (3.22) of Ref. [JTT 



Z^a 8 dr 
da T = 4vrC(3) — ) da s = yda T , (3.5) 

1 — T 

$ t (r, i/) = t~ ■ .„ F(l +iu,l- iv\ 2; z) 



1 + r) 3 sinh 7rz/ 



where 2; is defined in Eq. (|2.52|) . 
It gives 



= CJ z 1 ^ * 2 t (r, v) \h(A + mi , (3.6) 



r / , / " (1 x<1 j_ ^ dy ri j_ q 2n , 1 + 2 3 + 2x 2 

Mi)= /„ 1 + y -7^ — ^ = (1 + 3x ) In 2 - x 



y 7 2/(1 + 7/) 3 v 7 x 2 2(1 + x 2 ) 2 ' 

The properties of the function $ 2 (r, u), which determines the angular distribution of the 
produced ortho-Ps, are described in Ref. |Tl| . We mention here only that the typical values 
of r are r ~ 0.1 which correspond to a characteristic emission angle of ortho-Ps 

(3.7) 

The function fs(x) is large at small values of x 

f 3 (x) = 2 An - - at x < 1 (3.8) 

and small at large a; 

/3(a?) = at x > 1 • 

In the LLA we have 

h(x) + h(x) ^21n 7 2 , 2m e < £ <2m e7 . (3.9) 
For the spectrum of ortho-Ps in the LLA we obtain the expression 

ds r°° 

da ortho = 16C(3) z/V ln 7 2 — , B{v) = / $ 2 (r, 1/) dr . (3.10) 

£ JO 

Note that B(u) = 0.3770 for Au, 0.3654 for Pb and 

B[y) = 0.6137- 1.7073 v 2 at v < 1 . (3.11) 

After integrating the spectrum ( ft.lOD over £ in the region 2m e ~ e ~ 2m e 7 we obtain the 
total cross section in the LLA 
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Oortho = 8C(3) v 2 a B(v) In 7 



(3.12) 



The ortho-Ps production in lowest order in Za (production via three photons) is de- 
scribed by Eq. (|3.12|) with B(0) = 2 — 2 In 2 = 0.6137. Therefore, the relative order of the 
Coulomb correction (corresponding to the production with 5, 7, 9,. . . photons) is given by 
ratio 



B(0) 



(3.13) 



This ratio is equal to —39 % for the RHIC, for the LHC is reaches —40 % in the Pb-Pb 
mode and —3.6 % in the Ca-Ca mode. 



IV. DISCUSSION 

1. We have calculated the cross sections for the production of parapositronium (omitting 
terms ~ (Za) 2 /L 2 < 0.4%, L = hi7 2 ) and orthopositronium (omitting contributions ~ 
(Za) 2 /L < 4%) at the RHIC and the LHC. Our results are summarized in Table |T[ From 
that Table we observe that the production rate is well above 10 5 events per day, therefore 
measurements will be possible with a high precision at these colliders. We note that the 
production rate for the long living ortho-Ps is also relatively high, irrespectively of the 
suppression factor (Za) 2 . 

2. In these calculations we take into account not only the lowest order in the parameter 
Za but the whole series in Za — the so called Coulomb correction. They decrease the total 
cross section for the para-Ps production by about 30 % and 50 % for Au-Au and Pb-Pb 
collisions and by about 40 % for ortho-Ps production in the same collisions. For Ca-Ca 
collisions the CC are small ~ 3 %. 

3. We have presented the inclusive cross sections (see Eqs. ( j2.23| ), ( |2.28| ), (|2.62| ), ( |3.6|) 
and Figs. ||FJ) and the spectra for the Ps production (see Eqs. ( |2.45|) , (|2.56|) , (|3.10|) and 



Fig. |5|). From these results we conclude that the number of relativistic Ps with a Lorentz 
factors up to the Lorentz factor of the ion 7 will be large enough to be detected. 

4. The signature of the events is quite clear: the colliding ions practically do not change 
their motion and after the collision they remain in their bunches. In the central detector 
no hadrons should be observed. The further separation of the relativistic Ps from the back- 
ground can be performed using a method tested in a Serpukhov experiment (see Ref [0). 
In that experiment a special channel with a vacuum pipe of about 40 m length was used 
and a weak transverse magnetic field inside this pipe eliminates the charged particles. At 
the end of the pipe a strong magnetic field forces the Ps to break into e + and e~ which are 
easily detectable. 
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APPENDIX A: THE BORN AMPLITUDE FOR THE PARA-PS PRODUCTION 

It is useful to introduce the "almost light-like" 4-vectors 

M 2 M 2 ~ 9 M 6 ~ ~ 9 

Pi = Pi P 2 , P 2 = P 2 Pi, P? = ^r, s = 2P X P 2 w 2P X P 2 > M 2 

s s s z 

and to decompose the 4-vectors ft into components in the plane of the 4-vectors Pi and P 2 
and in the plane orthogonal to them 

6 i a 6 , 2 ftPi 2ftP 2 
^ = otiP 2 + A Pi + <2U > «t = , Pi = • (Al) 

s s 

The quantities a« and $ are the so called Sudakov variables. It is not difficult to estimate 
that 

N kjj ! N w _M_ < l^ 2 l <<; 1 (A2) 

fa 2 qi P 2 Am e M ' a 2 2q 2 P x Am e M ' V ; 

In the used c.m.s. the 4-vectors ftj_ have x and y components only 

ft_L = (0, q ix , Qjy, 0) = (0, qi_L, 0) , q 2 ± = -q- ± . 

The relativistic nuclei in the Ps production can be regarded as point-like and spinless 
particles. Therefore, the Born amplitude has the form (see Fig. [I] for the notations) 

where M^ v is the amplitude for the virtual process 7*7* — > Ps. Due to gauge invariance we 
have 

qt = (aiP a + P1P1 + qi±T = . 
This expression can be transformed to (using the estimates (|X2|)) 

(P x + P[YM^ = 2P?M IU/ « -^f± M M „ . 

Pi 

Analogously we find 

a 2 



which leads to the result 



M B = — ^f^- . A3 

Pia 2 q{q 2 
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Using the rules (2.11) from Ref. [11] for the transition from the pair production AA — > 
AAe + e~ to the process AA — > AA+ Ps, we obtain the Born amplitude for the production 
of the para-positronium in the quantum state n x Sq 

AUnZa) 2 m e o?l 2 
M B = -4 ^4 % x 



x-Tr 
4 



V- - qi + m e „ „ V- ~ & + m e „ 

qi± 7 T5 n <?2± + <?2_l_ 7 T5 k qi± 



(p. 



Qi 



(P- - <?2) 2 - ml 



where p + = p- = p/2. Taking the trace and using the estimates 

Am 2 



(p + mp s ) i7 5 
! we obtain 



h z a 7 s (qi ± x q 2 _L) ■ rii 



n 3 m P 



ni 



q\q\ 



|Pi| 



Am 2 — q\ — q 2 



2 ' 



(A4) 



From that Born amplitude we find the inclusive cross section for the production of para-Ps 
summed over quantum states n 



d 3 a B 

~¥p~ ~ 8(2tt 



1 oo 

~TE~2 X / \ M B\ 2 $(cLi± + <l2_L-p±)d 2 q 1± d 2 q2± 

T ) S n =l J 



(A5) 



Now we derive a useful expression for the virtualities q 2 . From the obvious relations 
(Pi - qi) 2 = M 2 or 2q 1 P 1 = sa x + M 2 ^ = q\ = sa 1 f3 1 + qf ± and 2q 2 P x = s(3 2 + M 2 a 2 = 
q\ = sa 2 f3 2 + q 2 ± we obtain 



*2 = ^ 



9 2 



q 2 ± - (M« 2 ) : 



1 - /3i 1 - «2 

Further, from the relation (gi + q 2 ) 2 = p 2 = Am 2 we get (using the inequalities (|A2|)) 

sa 2 /3 1 +p\ = 4m 2 e . 



(A6) 



(A7) 



If we introduce variables x = Mj3i/{2m e ) and x = Ma 2 /(2m e ) we immediately obtain from 
Eqs. 0-(^7D and inequalities Q the expressions (ggg)-(g]2§. 



APPENDIX B: CALCULATION OF J B FROM EQ. ( ^23| ) 



Let us introduce the two-dimensional vectors k, = qjj_/(2m e ) and the quantities V\ = 
k 2 + x 2 , V 2 = k 2 + x 2 , V3 — V1+V2 + I. Then the integral Jr from Eq. ( ^.23[ ) transforms to 



1 

7T 



(ki x k 2 ) 2 
{V x V 2 V^f 



S(k 1 + k 2 



Pi 

2m f 



d 2 k\d 2 k 2 . 



(Bl) 



Now we use the usual Feynman trick 

1 



(TO^s) 2 



5! 



<5(X^ «j — 1) IL oijdotj 
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and perform the integration over kj: 

Jb = 3r / dax 

Jo J o 



l ~ a i , ai^fl - a.i — a 2 )(2 — ai — a 2 f 
da 2 — 4 



(B2) 



D — (1 — cti)(l — 0:2)1" + (2 — ai — a 2 )[l — a± — a 2 + (1 — «2)x 2 + (1 — ai)x 2 ] 
where r is given in Eq. (gjp . Substituting 

2oi = (l + f)w, 2o 2 =(l-t)w 

we present Eq. ( |B2|) in the form 

3 z* 1 z* 1 (1 — t 2 ) 

J B = -r duu 3 {l-u)(2-u) 2 dV 



(B3) 



D = T 



1 -u + -u 2 (l -t 2 ) 
4 ; 



+ (2 - u) 



1 - u + x 2 + x 2 - -w(l - t)x 2 - + t)x 2 



which is equivalent to Eqs. ( p.28|) -( ^30|) . 
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TABLES 



Collider, mode 


Lorentz factor 


Luminosity 




7 = E/M 


C, cm 2 s _1 


RHIC, Au-Au 


108 


2 • 10 26 


LHC, Pb-Pb 


3000 


2 • 10 27 


LHC, Ca-Ca 


3750 


4 • 10 30 



TABLE I. Parameters of the discussed colliders. 



Collider, mode 


Opara 


para-Ps 


^ortho 


ortho-Ps 




mb 


events/sec 


mb 


events/sec 


RHIC, Au-Au 


8.7 


1.7 


6.8 


1.4 


LHC, Pb-Pb 


78 


155 


24 


48 


LHC, Ca-Ca 


0.41 


1650 


0.0086 


35 



TABLE II. Cross sections and event rates at the RHIC and LHC. 
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FIGURES 



Pj.E P/.Ej' 




FIG. 1. The amplitude M nn > for the production of Ps by n (n') virtual photons emitted by the 
first (second) nucleus. 
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FIG. 2. The inclusive cross section (|2.23[) as function of the scaled para-Ps energy e/(2m e 7) 
for the RHIC at different values of r = pj_/(2m e ) 2 . 
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FIG. 3. Same as Fig. | for the LHC (Pb-Pb mode) 
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FIG. 4. The electroproduction of Ps on a nucleus. 
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FIG. 5. The spectrum of the para-Ps as function of the scaled Ps energy e/(2m e 7). 
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